Disordered Optical Modes for Photon Management 
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Wave transport in disordered systems is a vast research topic, ranging from electrons in semicon- 
ductors, to light in random dielectrics, to cold atoms in laser speckles. In optics, light transport is 
conveyed by random electromagnetic modes and the wave can be localized about a point or extended 
throughout the system, depending on disorder strength, structural correlations and dimensionality 
of the system. Light localization phenomena are more dominantly present in two-dimensional sys- 
tems than in three-dimensional ones and their optical modes can be tailored to a greater extent. 
Here, we show that it is possible to make use of the properties of lower- dimensional disordered 
structures to obtain photon management in three-dimensional space. More particularly, we argue 
that two-dimensional disorder and wave interferences can be exploited to improve the performance 
of light absorbers or emitters. Our findings have direct applications for enhancing the absorption ef- 
ficiency of third-generation solar cells in a relatively simple and easy-to-implement way and provide 
unique opportunities for high extraction efficiency light-emitting diodes. 



In disordered optical materials, the multiple scattering 
of light and the interferences between propagating waves 
lead to the formation of electromagnetic modes with 
varying spatial extention and lifetime [lH3| . Their super- 
position gives rise to complex intensity speckle patterns 
that reflect the underlying transport process. Photon lo- 
calization 1 4], in analogy to the Anderson localization of 
electrons [5| , sets in when the scattering is strong or the 
dimensionality of the system is low. Two-dimensional 
disordered structures have been the playground over the 
years for the study of complex optical phenomena, in- 
cluding light localization [6-9] - two being the marginal 
dimension below which all states, in principle, are local- 
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ized [10] - and random lasing [11H13|. A greater con- 
trol over transport is achieved in two-dimensional sys- 
tems than in three-dimensional ones. Two-dimensional 
random structures in planar waveguides, in particular, 
constitute an important critical case that offers interest- 
ing possibilities. 

Already in the early publications on light localiza- 
tion [3, [HI , the absorption of light by a random medium 
was shown to be considerably enhanced by strong in- 
terference effects. Strongly scattering three-dimensional 
disordered structures appeared as a valuable candidate 
for designing strong light absorbers for photovoltaics ap- 
plications [16j . Nowadays, the main challenge taken up 
by the scientific community in the field of solar energy, 
however, consists in maximizing the absorption of solar 
radiation by absorbing layers much thinner than any pos- 
sible three-dimensional random structure. Thin-film so- 
lar cells can benefit a lot from the so-called photon man- 
agement, where a photonic structure is used to "trap" 
light in the absorbing medium [17]. Current strategies to 
improve their performance involve the use of plasmonic 
nanoparticles or gratings [l8l-[2o| , resonant dielectric par- 
ticles 21] or two-dimensional deterministic (periodic or 
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FIG. 1. Sketch view of a thin randomly-nanostructured film 
with thickness t, circular holes of diameter d and hole filling 
fraction /. Two-dimensional multiple scattering of light in 
the film leads to the formation of disordered optical modes to 
which coupling from free space is possible. 



quasi-periodic) photonic materials [22l-l25j. 



In this paper, we show that two-dimensional disordered 
patterns of air holes in planar dielectric waveguides pro- 
vide a conceptually new way to realize highly absorbing 
thin films. The structure under consideration is realistic 
and hence, the finite thickness of the film and the light 
coupling process between the in-plane optical modes and 
free space are intrinsically part of the physics involved. 
A delicate interplay between in-plane multiple scattering, 
out-of-plane losses, and optical absorption determines the 
optical properties of the system. For thin absorbing films, 
we find that it is possible to obtain a significant increase 
of the absorption over a very broad wavelength range, a 
wide angular range and for both polarizations of incident 
light - a crucial requisite for photovoltaic technologies. 
Although the coupling to this optical modes is mediated 
by interference effects, the highest absorption enhance- 
ment is surprisingly not obtained in the frequency range 
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in which they are the strongest, i.e. where the modes 
are spatially localized. The spectral and angular optical 
responses can be tuned by adding structural correlations 
in the disorder, which allows to match the wavelength 
and bandwidth of the absorption enhancement accord- 
ing to specific needs. The concept has the additional 
advantage that it is potentially simple to implement on a 
large scale and is compatible with current thin-film solar 
cell technologies. In addition to applications for light ab- 
sorbers, the same concept can be applied to enhance the 
extraction efficiency of light-emitting diodes for lighting 
applications, of great interest nowadays 26|, 27]. In that 
case, this allows to create an integrated optical source of 
diffusive light for which the angular profile can be tuned 
- a desireable property that is very difficult to obtain 
with conventional technology. 

The system under study is typically a thin dielectric 
membrane containing a random pattern of circular holes 
going through the entire film thickness (see Fig. [I]). The 
concepts presented below are valid for any type of two- 
dimensional random systems of finite thickness and can 
be applied to any inorganic or organic material. The 
choice of considering layers with small thicknesses is mo- 
tivated mainly by the growing demand in thin-film tech- 
nologies for photovoltaics and the need to enhance the 
light absorption for a higher photovoltaic efficiency at a 
lower cost. Let us first go briefly into the main physical 
concepts that describe light transport in such a struc- 
ture. Light is essentially confined to the plane of the 
film by index-guiding. Due to the presence of the holes, 
which serve as scattering centers, the light waves are 
multiply-scattered in the plane, giving rise a quasi- two- 
dimensional type of transport [28]. Because of the finite 
thickness of the film, however, the optical modes formed 
by the multiple scattering process are leaky and thus, 
allow for a coupling process with the third dimension 
to take place. The high density of states in the film 
makes such that the single scattering function is very 
much peaked into the plane, thereby leading to a very 
efficient light trapping even for planewaves at normal in- 
cidence on the slab. 

To demonstrate that this mechanism indeed leads to a 
very efficient broadband light absorber, we calculate the 
absorption spectrum for a planewave in air at normal in- 
cidence on an absorbing film patterned with circular holes 
at random positions, as in Fig.Q] The numerical analysis 
is performed by the three-dimensional Finite-Difference 
Time-Domain method [29] using periodic boundary con- 
ditions on a large area to mimic infinitely large struc- 
tures. For convenience, all dimensions in this study are 
given as a function of the film thickness t and can be 
rescaled according to one's need. Here, the circular holes 
have a diameter d = 2. It and filling fraction / = 30%. 
The material of the film has a permittivity e = 12 and 
an absorption length t{ = 33t, which, for the sake of the 
analysis, was considered constant over the wavelength 
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FIG. 2. Absorption spectra of a thin film including a random 
pattern of holes (blue line) and a bare thin film (black line). 
The inset shows the angular dependence of the absorption at 
t/X = 0.15 for incident light with s- and p-polarization (red 
circles and orange squares, respectively). The random pattern 
of holes leads to a large absorption of the incident light over 
a wide angular range and for both polarizations. 



range of interest. This corresponds approximately to the 
absorption of a film of amorphous silicon of thickness 
t = 75 nm at a free space wavelength A « 800 nm [30[ 
and compares well also with the absorption of some or- 
ganic materials used for thin-film solar cells. 

The absorption A of the patterned film is plotted in 
Fig. [2] together with the absorption of the bare film. We 
find that the random pattern of holes enhances signifi- 
cantly the absorption of the film over the entire wave- 
length range under consideration, the integrated absorp- 
tion being increased by a factor of about 5. This absorp- 
tion enhancement occurs in spite of the fact that 30% 
of absorbing material has been removed. The inset of 
Fig. [2] shows the angular dependence of the absorption 
at t/X = 0.15 for both s- and p-polarizations. One can 
see that the absorption remains particularly high even 
at large angles and for both polarizations of the incident 
light. This is a property of paramount importance for 
solar panels since it allows to retain a high efficiency also 
for diffuse illumination. 

To gain more insight into the mechanism behind this 
strong absorption enhancement, it is instructive to ob- 
serve how the electromagnetic energy density is dis- 
tributed in the nanostructured film. The repartition 
of the energy density between a dielectric material and 
free space allows to determine (at thermal equilibrium) 
the maximal enhancement of absorption in the mate- 
rial [I?], (HJ. In Fig. [3j the steady-state electromag- 
netic energy density distribution in the film is shown for 
two different frequencies (t/X = 0.09 and 0.15), consider- 



3 






Speckle intensity I 



FIG. 3. Top and side views of the electromagnetic energy den- 
sity inside a non-absorbing sample for an incoming planewave 
at frequencies t/X = 0.09 (left) and t/X = 0.15 (right). The 
normalization on the color scale is made according to the high- 
est value at the specific frequency. The polarization is along 
the x-direction and the incoming wavevector is at normal inci- 
dence to the plane of the slab, along the z-direction. The high 
energy density in the film implies an efficient light trapping, 
which explains the high absorption efficiency. 



FIG. 4. Speckle intensity distribution p(I) in the two- 
dimensional random medium at frequencies t/X = 0.09 
and 0.15. Strong light scattering and interferences between 
multiply-scattered waves give rise to slowly- decaying tails 
in the speckle intensity distribution that deviate from the 
Rayleigh law (black dashed line) . The variance of the distribu- 
tion (in the inset) shows that interferences between multiply- 
scattered waves are significant for transport over a broad 
wavelength range. 



ing again an incoming continuous planewave polarized 
along x at normal incidence on a non-absorbing pat- 
terned film with the same structural parameters as in 
Fig. [2j The remarkably high energy density in the film 
compared to free space (enhancement factors of about 
16 and 40 at tj X = 0.09 and 0.15, respectively) indicates 
a very efficient light trapping. The superposition of the 
excited in-plane modes is characterized by a spatial pat- 
tern dominated by strong fluctuations, inherent to the in- 
terferences between multiply-scattered waves. Different 
levels of spatial confinement can also be observed as a 
function of frequency. In leaky two-dimensional random 
media, the typical extention of the modes is expected to 
rely on the competition between wave localization effects 
and out-of-plane losses. 

The nature of the optical modes in the two-dimensional 
random medium can be further understood by analyzing 
the speckle distributions of the field and the intensity in- 
duced by the in-plane multiple light scattering. Their 
analysis allows to investigate how light propagates in the 
medium and how the coupling process with free space 
takes place. Since the transport is essentially in the 
plane of the film, this analysis can be made in purely two- 
dimensional systems. For the sake of the demonstration, 
we focus on the lowest-order TE-polarized propagating 
mode of the slab that is the most relevant for such thick- 
nesses. The refractive index of the dielectric material is 
set to the effective (frequency-dependent) refractive in- 



dex of the propagating mode as to account effectively for 
the finite thickness of the film. A large number of H z 
dipole sources with random phases are placed at random 
positions in the computational domain to excite all pos- 
sible modes. The spatial field maps recorded after long 
times then make it possible to retrieve important infor- 
mation on transport. 

The speckle intensity distribution, which describes the 
intensity fluctuations in random media, indicates the 
nature of the wave transport regime. In the case of 
an uncorrelated Gaussian disorder, the distribution of 
speckle intensity is expected to follow the Rayleigh law 
p(I) ~ exp(— I) with I = 1/(1) the speckle intensity nor- 
malized to its average value. Strong interference effects, 
such as light localization, typically lead to a strongly fluc- 
tuating intensity, characterized by a distribution with a 
slowly-decaying tail [3[. The speckle intensity in large 
purely two-dimensional structures, shown in Fig. |4] for 
t/X = 0.09 and 0.15, reveals deviations from the Rayleigh 
law at large intensities. The intensity fluctuations are 
larger at the lower of these two frequencies, that is in ac- 
cordance with the speckle patterns shown in Fig. [3l The 
variance of the distribution, shown in the inset of Fig. [4j 
indicates that strong interference effects are significant 
over a broad wavelength range. 

While the speckle intensity distribution gives infor- 
mation on the role of interferences on transport in the 
medium, a Fourier analysis of the fields makes it possi- 
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ble to infer the nature of the propagating modes in the 
random media and to extract the average dispersion of 
the system. The Fourier components of the electric field 
along the x- and ^-directions in the system are shown in 
Fig. [SK- The complex field pattern in real space is due 
to the superposition of the different optical modes in the 
medium. In reciprocal space, it is possible to separate all 
these modes, each point corresponding to a certain den- 
sity of (polarized) /c-space resonances [33]. For instance, 
the two lobes observed for both field components cor- 
respond to the most probable directions of propagation 
and their associated propagation constant. Such a study 
allows to investigate the nature of light coupling between 
different optical systems [34|. The coupling process be- 
tween in-plane modes and free space, of particular im- 
portance in this study, translates into a phase-matching 
condition that is met only for in-plane wavevectors k|| ly- 
ing within the so-called radiation zone, i.e. for |k|| | < /co, 
with ko = 2tt/X the wave vector in free space. 

By making use of the isotropy of the structure, we 
calculate the radial-averaged density of /c-space reso- 
nances M(t/\,k r ) = \E x (t/\k r )\ 2 + \E y (t/\k r )\ 2 with 
k r = |k|| | at different frequencies, see Fig. [SB, to describe 
the dispersion of the average propagating mode in the 
random medium. At lower frequencies, the distribution 
of resonances exhibits a well-defined peak, corresponding 
to a propagating mode with a certain effective refractive 
index. With increasing frequency, the wavelength in the 
medium becomes comparable to the size of the refractive 
index fluctuations due to the presence of the holes and, 
on the occurance of their broad Mie resonances, the peak 
broadens in reciprocal space [3|- In essence, the /c-space 
resonances lying in the radiation zone, delimited by the 
light line, are directly related to the in-plane light scat- 
tering process. This suggests that the coupling process 
can be controlled by tuning the structural parameters of 
the random medium. The /c-space resonances span the 
radiative zone quite uniformly, which explains why the 
coupling efficiency only depends weakly on the angle of 
incidence of light from the third dimension. Finally, a 
large part of the average mode lies out of the radiative 
zone, implying that light coupled to the film can remain 
within it for a certain time. In absorbing media, this 
trapping effect increases the probability for light to be 
absorbed. 

In summary, the distribution of the /c-space resonances 
of purely two-dimensional disordered structures is found 
to describe well how light propagates in the medium and 
how in-plane optical modes and free space are coupled to 
each other. This combination explains the strong absorp- 
tion enhancement that can be achieved by random pat- 
terning. As we will now show, it is possible, in addition, 
to tune the structural parameters of a random medium 
to optimize this absorption enhancement and calibrate 
its spectral response to a desired wavelength band. 

We have considered so far a random structure in which 



the position of the holes was uncorrelated (apart from 
a non-overlap constraint) and have seen that its opti- 
cal response does not exhibit sharp spectral and angular 
features. These optical characteristics, however, can be 
changed by adding correlations in the disorder. Struc- 
tural correlations, such as short-range order in amor- 
phous materials, provide an additional degree of freedom 
to control light transport and interferences in disordered 



materials [35|, [36|. In two-dimensional structures, corre- 
lations have been used to create large and isotropic pho- 
tonic band gaps [37[ as well as to gain a better control 
over random lasing processes [13]. To illustrate the pos- 
sibilities that structural correlations offer for light cou- 
pling purposes, we investigate the optical properties of 
an amorphous structure that possesses a short-range or- 
der in the position of the holes but lacks any long-range 
order (e.g. periodicity). 

We calculate the absorption spectrum of the two- 
dimensional amorphous structures for holes with the 
same diameter and at the same filling fraction, see 
Fig. [6]A. The effect of short-range correlations on the ab- 
sorption is striking. While the simple random structure 
provides a broadband absorption enhancement over the 
entire range of interest, the amorphous structure exhibits 
an absorption comparable to that of the bare film at lower 
frequencies, followed by a large absorption enhancement 
at frequencies t/X ~ 0.16. 

This behavior is perfectly explained by the distribu- 
tion of resonances shown in Fig. [6B for the amorphous 
case. At frequencies below to t/X « 0.07, a depletion of 
the density of /c-space resonances in the radiative zone is 
observed. Due to their very low number, light incident 
from free space on the film crosses it almost without be- 
ing scattered. This peculiar phenomenon is specific to 
our system of finite thickness and illustrates nicely that 
the coupling process is not a property of the single scat- 
terer but that of the disordered medium as a whole, in 
both the random and amorphous cases. 

This depletion of /c-space resonances is in fact accom- 
panied by two other remarkable effects, both of which are 
inherent to the short-range order in the structure and ev- 
idently do not require a long-range periodicity. First, at a 
frequency t/X « 0.07, corresponding to the Bragg condi- 
tion for which the wavelength in the medium A e = X/n e 
equals twice the typical distance a between nearest neigh- 
bours, X e = 2a, the dispersion undergoes a strong defor- 
mation. The lattice resonance is expected to enhance the 
in-plane scattering strength in this range [38]. Second, 
above this frequency, a branch with a negative slope is 
observed, strongly reminding the dispersion relations of 
periodically-nanostructured photonic crystals [39]. On 
the occurance of this new dispersion branch for k r = 0, 
the number of leaky resonances is strongly increased, 
thereby explaining the stronger absorption in this range. 
Furthermore, while the coupling to the random system is 
expected to be efficient at all angles, the folded branch 
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FIG. 5. (A) Distribution of the Fourier transform of the electric field components in the plane of the film, \E X \ and \E y \, at 
the frequency t/X = 0.15. Two-dimensional multiple scattering and wave interferences lead to the formation of complex optical 
modes, the superposition of which give complex speckle patterns. Light incident on such a random medium from the third 
dimension can couple to the structure through the /c-space resonances lying in the radiative zone (central disk delimited by a 
light blue solid circle). (B) Radial- averaged density of /c-space resonances M(£/A, k r ) in the two-dimensional random structure 
and the white lines are cross-cuts of the resonance distribution at frequencies t/A = 0.09, 0.12, 0.15 and 0.18. The light line is 
shown as a light blue dashed line. The leaky resonances in the radiation zone, above the light line, make the coupling between 
free space and in-plane modes possible. 



in the amorphous structure makes such that coupling 
at an angle may be more efficient in some frequency 
range. Disorder correlations therefore allow for a tun- 
ing of both the spectral and the angular light coupling 
of two-dimensional disordered media. For applications 
in lighting, one could take advantage of the negatively- 
sloped branch in amorphous structures to design light- 
emitting diodes with some iridescent diffuse light. All in 
all, structures with correlated disorder are promising can- 
didates for future integrated optical devices, as imposing 
only a weak constraint on the position of the scattering 
elements, such as a minimum distance, is enough to in- 
duce radically different optical characteristics. 

Finally, a critical parameter for the description of light 
absorption or extraction in thin films is the average life- 
time of the modes r. In particular, it is known that 
the absorption efficiency of nanostructured thin-film so- 
lar cells relies on a delicate interplay between the leak- 
age rate and the absorption rate of the material [17j. 
Typically, a good compromise in terms of efficiency and 
bandwidth is found when both rates are comparable. On 
the other hand, in light-emitting diodes, light should be 
extracted as fast as possible to avoid any possible re- 
absorption by the device 26, H?)- In this case, the av- 
erage lifetime of the modes should preferably be much 
smaller than the absorption time of the material. As 
the average lifetime of the modes in disordered systems 
depend intrinsically on the multiple scattering process, 
we expect it to be tunable by scattering strength and 



disorder correlations in thin two-dimensional disordered 
media. 

To investigate this possibility, we calculate the aver- 
age lifetime of the modes in the real (three-dimensional) 
random and amorphous structures by performing a tem- 
poral Fourier analysis of the decay of the field in the 
structure after excitation by many dipoles. In Fig. [7J the 
lifetime is compared with the total number of leaky res- 
onances in the system per unit frequency, evaluated as 
Ni (a/A) = J Jj k | <feo M(a/A, k||)dk|| from the resonance 
distribution of the corresponding two-dimensional struc- 
tures. At higher frequencies, the lifetime of the modes in 
the random and amorphous systems are very compara- 
ble, due to the fact that the number of /c-space resonances 
lying in the radiative zone is almost the same in both sys- 
tems. In contrast, at lower frequencies, the lifetime of the 
modes in the amorphous system becomes extremely long, 
as a result from the depletion by about an order of mag- 
nitude of the number of such resonances at frequencies 
close to t/X ~ 0.07. This phenomenon agrees with the 
previous analysis showing that the amorphous material in 
this frequency range can be transparent to light incident 
from the third dimension. In essence, this shows that 
the average lifetime of the modes in the system can be 
controlled by its structural parameters, thereby opening 
the way for the design of innovative disordered photonic 
structures with specific optical properties. 

In conclusion, controlling the light coupling process be- 
tween free space and thin or ultrathin films is a topic of 
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FIG. 6. (A) Absorption spectrum of a two-dimensional amor- 
phous structure (red solid line), compared to the absorption 
spectrum of bare film (black line) and of the random struc- 
ture (blue dashed line). The modelled amorphous structure 
is shown in the inset. (B) Radial- averaged density of /c-space 
resonances M(£/A,&v), as in Fig. [4] in the two-dimensional 
amorphous structure. The amorphous structure, due to struc- 
tural correlations, leads to a larger increase of the absorption 
on a smaller spectral range compared to the simple random 
case. The reduced and enhanced absorption are explained by 
the depletion and increase of the density of resonances for a 
planewave at normal incidence on the film at lower and higher 
frequencies, respectively. 



great interest for many applications, in particular in the 
field of energy efficiency. In this paper, we have proposed 
an original way to take up this challenge, namely that of 
using the optical modes formed by the multiple scatter- 
ing of light in two-dimensional disordered media of finite 
thickness. We have shown in particular that the strength 
of the in-plane scattering and the structural correlations 
in the disorder offer the possibility to tune the spectral 
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FIG. 7. Average lifetime of the modes in the two-dimensional 
random (blue dashed line) and amorphous (red solid line) 
structures. The inset shows the number of leaky resonances 
per unit frequency in the random and amorphous structures. 
The depletion in the number of leaky resonances in the amor- 
phous structures at frequencies t/X « 0.07 leads to a remark- 
able increase of the average lifetime of the modes. 

and angular optical properties of the light coupling pro- 
cess, as well as the average lifetime of the modes. Such 
disordered patterns are likely to be easier to implement 
on an industrial scale compared to their periodic coun- 
terpart, as no long-range order is required and the optical 
properties of the material are far more robust to struc- 
tural imperfections. For these reasons, we believe that 
two-dimensional disordered materials could form a new 
generation of high-efficiency (inorganic or organic) thin- 
film photovoltaic devices and light sources. 

By bringing the physics of random media into the field 
of energy efficiency, this study also opens the route to- 
wards the use of complex optical phenomena of interest 
so far for fundamental science mainly, such as light lo- 
calization, for future technologies. For instance, light 
emitters (e.g. quantum dots) embedded in such two- 
dimensional disordered media of finite thickness would 
see their rate of spontaneous emission dramatically in- 
crease due to a very strong light-matter interaction, as 
observed recently in disordered one-dimensional photonic 
crystal waveguides [40]. By a clever disorder engineering, 
one may possibly boost up the light emission efficiency 
in thin films in specific frequency ranges. 

Finally, let us emphasize the fact that the concepts 
presented in this article rely on the properties of waves 
in general and thus, are not restricted to electromagnetic 
waves in the optical range. Thin two-dimensional dis- 
ordered structures could be designed to operate for mi- 
crowaves or sound waves, and be used for shielding appli- 
cations. Two-dimensional random media of finite thick- 
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ness are rich of new physical concepts for wave manage- 
ment in three-dimensional space. This article paves the 
way for more research in this direction as well as for the 
development of concrete applications. 

Methods 

The random structures, on one hand, were gener- 
ated by Random Sequential Addition [4l| of hard cir- 
cular cylinders at a filling fraction of 30% on a square 
area with periodic boundary conditions. The amorphous 
structures, on the other hand, were generated by the 
Lubachevsky-Stillinger algorithm on the same area [42|. 
This algorithm was used to impose short-range disorder 
correlations in the position of the cylinders. To avoid the 
formation of large ordered clusters, as in polycrystalline 
structures, a polydispersity of 1.1 in the size of the cylin- 
ders was used. The growth of the cylinders was stopped 
at a filling fraction of 69% and their diameter was set 
to a constant diameter such that the final filling fraction 
would be 30%. 

For the refractive index of the dielectric medium in 
the two-dimensional calculations, we used the effective 
index of refraction of the lowest-order TE propagating 
mode, calculated by an online mode solver for multi- 
layered media [43| at different wavelengths. A polyno- 
mial fit was used to model the refractive index of the 
dielectric film over the entire range of interest: n s = 
3.6446 - 0.1074 (f ) - 2.213 10" 3 (f f + 1.744 10~ 4 (f ) 3 . 
For the effective refractive index of the two-dimensional 
randomly patterned array, the two-dimensional Maxwell- 

Garnett mixing rule can be used: n e = n s ^Jl + ilf a •> 

2 _ 2 

with a = n %~, n % with n c = 1 the refractive index of 
the cylinders and n s the effective refractive index of the 
lowest-order TE propagating mode defined above. 

Simulations were performed with the Finite-Difference 
Time-Domain method [29| using a freely available soft- 
ware package [lif. The mesh resolution was set to have 
at least 10 grid points per wavelength in the higher re- 
fractive index material. Periodic boundary conditions in 
the lateral directions (x and y) were used in all simu- 
lations. The lateral size of the computational domain 
was L = 200t in the two-dimensional calculations and 
L = 50t in the three-dimensional calculations. The ab- 
sorption A was calculated by recording the reflected and 
transmitted fluxes, R and T, respectively, for long times 
under excitation by planewave pulses at different frequen- 
cies, and calculating A = 1 — T — R. For the studies 
on the two-dimensional speckle field and intensity distri- 
bution of the random and amorphous structures, a to- 
tal of 2000 dipoles, uniformly distributed throughout the 
entire computational domain and with random phases, 
was used. For the lifetime calculations in the three- 
dimensional structures, 640 E x , E y and E z dipoles were 



placed at random positions in the film and the decay of 
the electric field was recorded at 20 different points uni- 
formly distributed in the film. 
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